Testing the occurrence of Late Jurassic true polar wander using the La Negra volcanics of northern Chile

True polar wander (TPW) is the reorientation of the crust-mantle system driven by the redistribution of masses in the mantle and on the Earth’s surface. In the ideal case, characterization of TPW requires paleomagnetic constraints on the motion of all major plates and independent reconstructions of relative plate positions. While such complete datasets are absent for pre-Mesozoic TPW inferences due to the absence of oceanic plates, they are available for the Late Jurassic (165-145 Ma) “monster shift”, a ∼30◦ amplitude proposed TPW event. Here we perform paleomagnetic sampling and Ar-Ar geochronology on the La Negra volcanics of Northern Chile, producing two new paleomagnetic poles with ages 165.8 ± 1.8 Ma (1σ; 84.3◦N 0.9◦E; α95 = 7.6◦; N = 28) and 152.8 ± 0.8 Ma (84.5◦N 256.4◦E; α95 = 10.8◦; N = 18). By combining these data with other recently published results, we compute a net lithospheric rotation of 25.3◦ ± 7.3◦ (1σ ) at a mean rate of 1.21◦ ± 0.35◦ My−1 between 170 and 145 Ma with a peak rate of 1.46◦ ± 0.65◦ My−1 between 160 and 145 Ma. These rates are consistent with inferences from the Pacific Plate, implying true whole lithosphere rotation. Given coherent motion involving the entire lithosphere, we conclude that the Earth underwent rapid TPW between approximately 165 and 145 Ma, potentially driven by the cessation of subduction along the western North American margi

Revisiting tectonic corrections applied to Pleistocene sea-level highstands

Tectonic displacement contaminates estimates of peak eustatic sea level (and, equivalently, minimum continental ice volumes) determined from the elevation of Quaternary interglacial highstand markers. For sites at which a stratigraphic or geomorphic marker of peak Marine Isotope Stage (MIS) 5e sea level exists, the standard approach for estimating local tectonic uplift (or subsidence) rates takes the difference between the elevation of the local highstand marker and a reference MIS 5e eustatic value, commonly chosen as +6 m, and divides by the age of the marker. The resulting rate is then applied to correct the elevation of all other local observed sea-level markers for tectonic displacement, including peak highstands of different ages (e.g., MIS 5a, MIS 5c and MIS 11), under the assumption that the tectonic rate remained constant over those periods. This approach introduces two potentially significant errors. First, the peak eustatic value adopted for MIS 5e in most previous studies (i.e., +6 m) is likely incorrect. Second, local peak sea level during MIS 5e is characterized by significant departures from eustasy due to glacial isostatic adjustment in response to both successive glacial–interglacial cycles and excess polar ice-sheet melt relative to present day values. We use numerical models of glacial isostatic adjustment that incorporate both of these effects to quantify the plausible range of the combined error and show that, even at sites far from melting ice sheets, local peak sea level during MIS 5e may depart from eustasy by 2–4 m, or more. We also demonstrate that the associated error in the estimated tectonic rates can significantly alter previous estimates of peak eustatic sea level during Quaternary highstands, notably those associated with earlier interglacials (e.g., MIS 11)

Rapid postglacial rebound amplifies global sea level rise following West Antarctic Ice Sheet collapse

Earth_Model_Data is a zipped folder containing the Earth model data for both the standard model (V3D_SD and V3D_RH). A readme file is in this folder. FPRINT_CODE is a zipped folder containing the fingerprint code. A readme file for the code is also in this folder. WAmask_512.gz is a mask for West/East Antarctica, used for masking out changes in East Antarctica. All other files are sea-level outputs for each of the runs on a degree 512 Gauss-Legendre grid (uniform longitudes but unevenly spaced latitudes, as described in the readme for FPRINT_CODE). Files are named SLt_??? and numbered from 0 (elastic response) to 105 (10 ky). They have a 1D layout, with the first line being the time tag in years followed by 512*1024 row entries. A time array is included (tt_v10.dat). More details are in README.tx

Pyrite-walled tube structures in a Mesoproterozoic sediment-hosted metal sulfide deposit

Unusual decimeter-scale structures occur in the sediment-hosted Black Butte Copper Mine Project deposit within lower Mesoproterozoic strata of the Belt Supergroup, Montana. These low domal and stratiform lenses are made up of millimeter-scale, hollow or mineral-filled tubes bounded by pyrite walls. X-ray micro−computed tomography (micro-CT) shows that the tube structures are similar to the porous fabric of modern diffuse hydrothermal vents, and they do not resemble textures associated with the mineralization of known microbial communities. We determined the sulfur isotopic composition of sulfide minerals with in situ secondary ion mass spectrometry (SIMS) and of texture-specific sulfate phases with multicollector−inductively coupled plasma−mass spectrometry (MC-ICP-MS). The sedimentological setting, ore paragenesis, sulfur isotope systematics, and porosity structure of these porous precipitates constrain the site of their formation to above the sediment-water interface where metalliferous hydrothermal fluids vented into the overlying water column. These data constrain the geochemistry of the Mesoproterozoic sediment-water interface and the site of deposition for copper-cobalt-silver mineralization. Metals in the hydrothermal fluids titrated sulfide in seawater to create tortuous fluid-flow conduits. Pyrite that precipitated at the vent sites exhibits large sulfur isotope fractionation (>50‰), which indicates a close association between the vents and sulfate-reducing microbiota. In the subsurface, base metal sulfides precipitated from sulfide formed during the reduction of early diagenetic barite, also ultimately derived from seawater. This model suggests dynamic bottom-water redox conditions at the vent site driven by the interplay between sulfate-reducing organisms and metalliferous fluid effluence

Snowball Earth climate dynamics and Cryogenian geology-geobiology

Geological evidence indicates that grounded ice sheets reached sea level at all latitudes during two long-lived Cryogenian (58 and ≥5 My) glaciations. Combined uranium-lead and rhenium-osmium dating suggests that the older (Sturtian) glacial onset and both terminations were globally synchronous. Geochemical data imply that CO2 was 102 PAL (present atmospheric level) at the younger termination, consistent with a global ice cover. Sturtian glaciation followed breakup of a tropical supercontinent, and its onset coincided with the equatorial emplacement of a large igneous province. Modeling shows that the small thermal inertia of a globally frozen surface reverses the annual mean tropical atmospheric circulation, producing an equatorial desert and net snow and frost accumulation elsewhere. Oceanic ice thickens, forming a sea glacier that flows gravitationally toward the equator, sustained by the hydrologic cycle and by basal freezing and melting. Tropical ice sheets flow faster as CO2 rises but lose mass and become sensitive to orbital changes. Equatorial dust accumulation engenders supraglacial oligotrophic meltwater ecosystems, favorable for cyanobacteria and certain eukaryotes. Meltwater flushing through cracks enables organic burial and submarine deposition of airborne volcanic ash. The subglacial ocean is turbulent and well mixed, in response to geothermal heating and heat loss through the ice cover, increasing with latitude. Terminal carbonate deposits, unique to Cryogenian glaciations, are products of intense weathering and ocean stratification. Whole-ocean warming and collapsing peripheral bulges allow marine coastal flooding to continue long after ice-sheet disappearance. The evolutionary legacy of Snowball Earth is perceptible in fossils and living organisms

Geobiology of a lower Cambrian carbonate platform, Pedroche Formation, Ossa Morena Zone, Spain

The Cambrian Pedroche Formation comprises a mixed siliciclastic-carbonate succession recording subtidal deposition on a marine platform. Carbonate carbon isotope chemostratigraphy confirms previous biostratigraphic assignment of the Pedroche Formation to the Atdabanian regional stage of Siberia, correlative to Cambrian Series 2. At the outcrop scale, thrombolitic facies comprise ~. 60% of carbonate-normalized stratigraphy and coated-grains another ~. 10%. Petrographic point counts reveal that skeletons contribute at most 20% to thrombolitic inter-reef and reef-flank lithologies; on average, archaeocyath clasts make up 68% of skeletal materials. In contrast, petrographic point counts show that skeletons comprise a negligible volume of biohermal and biostromal thrombolite, associated nodular carbonate facies, and ooid, oncoid and peloid grainstone facies. As such, archaeocyathan reefal bioconstructions represent a specific and limited locus of skeletal carbonate production and deposition. Consistent with data from coeval, globally dispersed lower Cambrian successions, our analysis of the Pedroche Formation supports the view that lower Cambrian carbonates have more in common with earlier, Neoproterozoic deposits than with younger carbonates dominated by skeletal production and accumulation. © 2013 Elsevier B.V.Jessica R. Creveling, David Fernández-Remolar, Marta Rodríguez-Martínez, Silvia Menéndez, Kristin D. Bergmann, Benjamin C. Gill, John Abelson, Ricardo Amils, Bethany L. Ehlmann, Diego C. García-Bellido, John P. Grotzinger, Christian Hallmann, Kathryn M. Stack, Andrew H. Knol

The sea-level fingerprint of a Snowball Earth deglaciation

Cap dolostones are thought to represent deposition from seas transgressing over formerly glaciated continental margins during Marinoan Snowball deglaciation. Nevertheless, facies associations within some cap dolostones indicate that an episode of regional regression punctuated these transgressive sequence tracts. To date, inferences of sea-level change during and after the Marinoan Snowball deglaciation have been interpreted using simple, qualitative arguments. In the present study, we explore the full spatio-temporal variability of sea-level change during Snowball deglaciation and its aftermath using a gravitationally self-consistent theory that accounts for the deformational, gravitational and rotational perturbations to sea level on a viscoelastic Earth model. The theory is applied to model Marinoan Snowball deglaciation on a generalized Ediacaran paleogeography with a synthetic continental ice-sheet distribution. We find that sea-level change following a synchronous, rapid (2 kyr) collapse of Snowball ice cover will exhibit significant geographic variability, including site-specific histories that are characterized by syn-deglacial sea-level fall or stillstand. Moreover, some sites that experience syn-deglacial transgression will continue to experience transgression in the post-deglacial phase. Taken together, these results suggest that sea-level change recorded by strata capping Snowball glaciogenic units may reflect a more complicated trajectory than previously thought, including deposition that was not limited to the deglaciation phase. These simulations, as well as others that explore the response to asynchronous melting and deglaciation phases of longer duration, demonstrate that an episode of regional regression interrupting a cap dolostone transgressive sequence tract may reflect one of several processes (or their combination): (1) near field adjustment associated with rapid local melting during an otherwise global hiatus in deglaciation; (2) post-glacial uplift of sites during a period of slowing deglaciation, and (3) a transition, at some sites, from a sea-level fall dominated by post-glacial uplift to a phase of sea-level rise due to eustasy and peripheral bulge subsidence throughout an extended (order 50 kyr or greater) Snowball deglaciation

Cap carbonate platform facies model, Noonday Formation, SE California

The Neoproterozoic outcrop belt of the Death Valley region, California, preserves an oblique cross section of the Noonday Formation, a mixed carbonate-siliciclastic platform that hosts distinctive basal Ediacaran cap carbonate−affiliated sedimentary structures, stromatolite textures, and δ^(13)C_(carb) values. The Noonday platform encompasses two depositional sequences that reveal two cycles of relative sea-level change within strata conventionally considered to record a single, rapid, postglacial sea-level rise. In updip localities, facies of the first depositional sequence record the transition from a carbonate ramp to a stromatolite-bearing, “tubestone”-textured, reef-rimmed platform; downdip, localities seaward of the reefal escarpment variably preserve a thin and condensed onlapping foreslope wedge. Base-level fall exposed the reef crest to karstic dissolution and propagated submarine incised valleys into the seaward margin of the reef. Overlying strata record the backfilling of a submarine incised valley and reestablishment of a back-stepping, carbonate-dominated ramp prior to a second subaerial exposure event that defines the contact between the Noonday and Johnnie formations. We address the relative contributions of syndepositional tectonism and recovery from low-latitude deglaciation in dictating Noonday platform architecture and the intra−Noonday Formation sequence boundary. Noonday Formation deposition coincided with extension of the Laurentian margin during disaggregation of the Rodinian supercontinent. Within this framework, previous work has suggested that the intra−Noonday Formation sequence boundary records growth faulting that reinforced differential topography, uplifting reef-rimmed horsts—exposing the reef crest to karstic dissolution—and downdropping grabens. However, we trace the intra−Noonday Formation sequence boundary seaward of the reef crest and demonstrate that, for a time, wave base was situated downdip of the reef escarpment on putatively downdropped fault blocks. Thus, if the Noonday margin were undergoing extension, then the creation of the intra−Noonday Formation sequence boundary required a concomitant decrease in accommodation due, perhaps, to postglacial isostatic uplift attendant with low-latitude deglaciation. We speculate that Noonday Formation sequence architecture records (1) immediate deglacial flooding, (2) shoaling and exposure due to isostatic rebound induced by either a hiatus in meltwater flux or rapid ice-sheet collapse against a background of global deglaciation, and (3) resumed flooding following complete deglaciation. As rift-related tectonism could amplify or counter glacial isostasy, inferences of the amplitude of local postglacial sea-level change will require robust estimates of syndepositional extension across the Noonday margin

Dynamic Time Warping of Paleomagnetic Secular Variation Data

We present and make publicly available a dynamic programming algorithm to simultaneously align the inclination and declination vector directions of sedimentary paleomagnetic secular variation data.This algorithmgenerates a library of possible alignments through the systematic variation of assumptions about the relative accumulation rate and shared temporal overlap of two or more time-series. The paleomagnetist can then evaluate this library of reproducible and objective alignments using available geologic constraints, statistical methods, and expert knowledge.We apply the algorithm to align previously (visually) correlated medium to high accumulation rate northern North Atlantic Holocene deposits (101– 102 cm/ka) with strong radiocarbon control. The algorithm generates plausible alignments that largely conform with radiocarbon and magnetic acquisition process uncertainty. These alignments illustrate the strengths and limitations of this numerical approach